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ABSTRACT 

In this report, the basic heat transfer and isothermal flow friction 
characteristics for five perforated nickel matrices, one perforated nickel 
fin with solid splitter matrix, and one solid nickel matrix, all of 
"parallel plate'’' configuration are presented. In addition, similar in- 
formation is provided for a 20° skew matrix fabricated from perforated 
nickel. Finally, performance of two offset rectangular fin matrices 
constructed from stainless steel, is evaluated. Where meaningful, a 
comparison of the performance of these surfaces is made. All matrices 
are investigated by the transient method of testing utilizing the maximum 
slope analysis of Locke and include the effects of longitudinal thermal 


conduction as reported by Howard. 
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l. INTRODUCTION 

While it is appreciated that the designer of a compact heat exchanger 
is confronted with many qualitative evaluations, nevertheless he is general- 
ly working specifically within limitations relating to flow frontal area, 
weight, pressure drop and volume. It is for this reason that different 
heat transfer surfaces may be effectively compared on graphs using co- 
ordinates of f and j (or f/j ratio) versus Np: 

As in many other disciplines, model theory is the most satisfactory 
approach from a cost effectiveness standpoint. The use of such methods 
as steady state "steam to air" tests or "transient" testing techniques 
provides the basic flow friction and heat transfer characteristics without 
necessitating the construction of a full scale heat exchanger. 

The data presented herein was obtained by the use of a transient test 
technique or "single blow". £7] This technique is one of heating a sample 
core to some predetermined value and then after thermal equilibrium is 
achieved in the specimen, instantaneously removing the upstream heat source 
and monitoring the downstream temperature response of the fluid at the core 
exit. 

The objectives of this thesis specifically were: to modify an experi- 
mental test rig to operate with a larger capacity blower and to test sever- 
al offset rectangular fin, of stainless steel, and "modified parallel" plate, 
of perforated nickel, sample cores. 

In the testing of perforated nickel specimens, the effect of varying 


shape and size of perforations was to be investigated. 





ae SUMMARY OF THEORY 

The “starting stage of operation", "single blow", or "transient test 
technique", as it has been called by various authors, had its inception 
with the analytical treatment by Anzelius in 1926. Nusselt in 1927, Hausen 
in 1927 and 1929, and finally Schumann in 1929 treated the problem of ana- 
lytically determining the fluid and solid temperature as a function of 
position and time for a porous solid, initially at a uniform temperature, 
when subjected to a change in the entering fluid temperature. Some repeti- 
tion occurred as the authors were unaware of the work performed by their 
contemporaries at the time [9] . Anzelius and Schumann made similar as- 
sumptions; however, the results are different. Anzelius simplified the 
problem; whereas, Schumann proceeded to obtain an explicit mathematical 
solution for the temperature distribution. This solution was for a liquid 
flowing through a porous solid and consisted of two infinite series which 
were related to Bessel functions. Schumann states [17] that if the fluid 
is a gas, the solution to the problem is more formidable. He proceeds to 
relate how dimensional analysis may be used so that his investigation re- 
mains usable for compressible fluids. 

The porous solid is considered to be homogeneous and at a uniform 
temperature throughout with the fluid flowing at this same temperature. 
With a sudden change in the entering fluid temperature, the temperature 
of the fluid and the solid as a function of both position and time re- 
mains to be evaluated. The theoretical generalized heating curves are 
developed from an energy balance on the porous solid as presented in [19] 
and restated in [1]. Inasmuch as experimental procedure used in this 
report is dependent upon the assumptions and boundary conditions of this 
energy balance, it is briefly summarized below: 
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Mept. > hb(tp-ts)dx = they (+, + BF ax) 


St > KET > a (885) 


b~—— dx ——m} 


Assumptions made in the analysis are: 


a. Properties of the fluid are temperature independent 
b. Fluid flow is steady 
c. Porous solid is homogeneous ; 


d. Thermal conductivity of both fluid and solid is infinite 
perpendicular to the fluid flow direction 


e. Thermal conductivity of fluid is zero in the flow direction 


Boundary conditions are [13]: 


Heat 
Heat 
Heat 
Heat 


Thus 


(a) The matrix is initially at a uniform temperature 

(b) At time equal zero, the temperature of the entering fluid 
changes instantaneously to a different, constant value, 
i.e., a step change in fluid temperature 

(c) The matrix boundaries are adiabatic 


2 
absorbed by solid = ACs ots dx 
MA Oss BY) 


transferred to the solid by convection = hb (t,-t.)d 


transferred from the fluid by convection = ac dx 
x 
transferredfrom solid by conduction = -k ss § as dx 


an energy balance for the fluid and solid o-. 
(1) eg Ascs os dx = ksAs Sts dx + hb (tp-t¢) dx 


(2) mop St dx + hb (tp-t,) dx =O 





Rearranging the above equations and substituting in the following para- 


meters {3] : 


W, 
T = generalized time variable = hA je- ae x | 
s Ss ee 


where h = unit conductance for convective heat transfer, 
Btu/hr ft °F 
A = matrix heat transfer area, £t2 
W.c. = matrix capacity, Btu/ °F 
= time, hrs. 
W. = mass of fluid contained within the matrix, lb. 
m = mass flow rate of fluid, lb/hr. 


x = distance from entrance of matrix along flow passage 
in the direction of flow, ft. 


L 


matrix flow length, ft. 


Regrouping the above equation for 4 and introducing the specific heat of 


the fluid, oe we get: 7 = HA @) = hA We Cp x 
Ws Cg m Ws Cs Cp L 


Since the thermal capacity of the fluid contained at any time within the 
matrix (Wee) is much less than the thermal capacity of the matrix (Woe), 


the second term can be neglected; accordingly: 








Te hA 6 (3) 
Ws Cs 

z = generalized position variable =hA x 

Mop L 
Inasmuch as the number of transfer units (Nou) = hA » this simply 

m Cp 
becomes 
s x 
z= Neu arias (4) 





A = kA (5) 





8 
me L 
P 
where io = matrix thermal conductivity, Btu/hr yee a Pret. 
A = solid matrix cross sectional area available for thermal conduc- 


$ tion, ft. 


m, Cs and L are as given above. 


Equations (1) and (2), together with (3), (4), and (5) yield: 


dts = yt 

ao le yor = Cre, (6) 
See 

a ap ts ike (7) 
Ze 


where t. = temperature of the solid 
t- = temperature of the fluid 
For the case of thermal conduction in the solid parallel to the direction 


of flow equal to zero, equations (6) and (7) simplify: 


ones - 
= = t, te (8) 
i SE a. 3) (9) 


It is with these two differential equations, that Schumann utilizing 
Bessel functions determined the time - temperature history for the solid 
after a step change in fluid temperature. For a more thorough treatise, 
the reader is referred to references [13] or [17] . 

The task of determining the ae value for a test point required the 


superposition of the experimentally obtained time - temperature heating 





Or cooling curve on the family of theoretical curves. Although this would 


uniquely determine the N, value, it can be appreciated as being a tedious 


Tu 
effort. 

Fortunately, in 1950, Locke [13] developed a "maximum slope" technique 
which bypasses this task. He differentiated the theoretical curves and 
plotted the slope of the curve versus v My, with Nou as a parameter. The 


expression for the slope is 
Fee t: Ze. 
deel = Nru_ -J, (2:1 Nru a Nat) 
dT /nru] { Nr 


where the downstream fluid temperature is evaluated at x = L so that z = 
Neu and te = Ceo: 

In the efforts of both Locke and Schumann, the thermal conductivity 
of the solid parallel to the flow was taken as zero. Although the effect 
of longitudinal conduction is not pronounced at high Reynolds number range, 
it is considerable at low laminar Reynolds numbers. Several authors have 
considered the effects of longitudinal conduction, introducing the conduc- 
tion parameter, A ° When AF © , the maximum slope is not a unique 
function of Neu but depends on . as well. Howard [7] has tabulated maxi- 
mum slopes as functions of Nou and /\ » as well as presenting graphs, 
which he evaluated by a finite difference technique with the aid of a 
digital computer. The evaluation of Neu in this report is based on these 
tabulated values in so far as practicable. 

Satisfying the aforementioned idealizations, initial and boundary 
conditions, leaves the determination of the conduction parameter, Ne , 
and the experimental evaluation of the maximum slope of the heating or 


cooling curve as the only two remaining items necessary to determine the 


heat transfer characteristics of the matrix. 





3. EXPERIMENTAL TECHNIQUES 


The results of Howard's finite difference calculations in reference 


[7] were used in this report for evaluating _ and ultimately the 
u 


Colburn j factor. For this reason, the idealizations imposed by Howard 
were restrictions placed on the experimental set up used. Howard L7] 
considered transient heating of a porous solid when subjected to a step 
change in fluid temperature including the material longitudinal conduc- 
tion effects. The idealizations are restated below: 


"(1) The flow entering the matrix is steady and uniform in 
velocity and temperature, and remains steady and uniform 
at any cross-section as it flows through the matrix. The 
thermal conductivity of the matrix is infinite in the 
direction normal to the flow and finite in the direction 
parallel to the flow. Thus the problem is one-dimensional 
in space. 


(2) The thermal capacity of the fluid contained at anytime 
within the matrix is small compared with the thermal 
capacity of the matrix. That is, the fluid is normally 
restricted to a gas and the results to be presented 


-will not be expected to apply for liquid fluids. This 
also means that for the fluid there will be no time de- 


pendent terms in the equations. 


(3) The thermal properties of the fluid and matrix are con- 
stant and uniform. 


(4) The convective heat transfer coefficient is some suitable 
average and remains uniform and constant. 


(5) At time zero the change in fluid temperature will be a 
step function, with the matrix and its entrained fluid 
initially at some uniform and equal temperature.” 
Idealization (1) was achieved, as far as steadiness and uniformity of 
temperature and velocity were concerned, by responsible experimentation. 
In an effort to improve on the previously used test apparatus of Bannon 


[1] , @ new inlet cone was designed and the test section relocated. The 


inlet cone was fashioned according to the specifications of Smith and 





Wang [18] for contracting cones giving uniform velocity at the throat. 
Although this design was specifically for circular cones, the epediileace 
square inlet cone fabricated, together with wire screen type flow straight- 
eners provided the required uniform velocity profile. Velocity profiles, 
taken at the matrix inlet over the entire flow range, verified required 
uniformity. The test section was placed upstream of the air metering 
device so that full advantage of the uniform velocity profile could be 
realized. Uniform temperature profiles within + 0.5 °F were obtained by a 
wire screen type heater explained in detail in APPENDIX B. The latter part 
of the first idealization pertaining to the thermal conductivity of the 
matrix cannot be achieved precisely; however, the effects of deviating 

from it can be made negligible in the solid, if resistance to heat transfer 
is small within the solid compared to the heat transfer resistance be- 
tween solid and gas. In other words, "transverse conduction" effects have 
been neglected; whereas, the effects of longitudinal conduction are consid- 
ered with the use of A. , the conduction parameter. In a matrix where 
adequate mixing takes place, no criterion is necessary for "transverse 
conduction" effects. 

The second idealization is self explanatory, i.e., the fluid is re- 
stricted to a gas. 

Idealization (3) for the thermal properties can be guaranteed if the 
temperature differences are kept small. A maximum air temperature varia- 
tion from 70 °F to 90 °F results in a maximum deviation from the arith- 
metic average in the order of + 1.5% in viscosity and is negligible in 
specific heat. 

The convective heat transfer coefficient presented as the final result 
of this experiment is an overall quantity; accordingly, idealization (4) 
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is satisfied. 

Satisfying the fifth idealization of a step change in fluid tempera- 
ture is perhaps the most rigorous requirement of all. Several investiga- 
tors have attempted to meet this requirement by various schemes. The 
technique of Bannon [1] and Theoclitus and Eckrich (20] appears to be 
most satisfactory; therefore, it was the one selected for this investiga- 
tion. 

Figure 1 displays a schematic of the experimental set up. Figure 2, 
is a photograph of the complete test apparatus. The system is comprised 
of an air source, flow straightener screen section, heater section, matrix 
test section, air metering device, and pressure and temperature measuring 
Systems. Each of these are discussed at length in APPENDIX B. 


Air was drawn through the "wind tunnel" by a Spencer Turboblower. 


The 0.0031" diameter electric wire heaters upstream of the matrix were 


energized so as to heat the air to roughly 20 degrees above ambient. The 
upstream and downstream fluid temperatures at the matrix test section were 
monitored to ascertain when the matrix had been heated to a uniform tempera- 
ture. When this condition was satisfied, the heaters were quickly shut 
down and the time-temperature history of the entering and exit air at the 
matrix was recorded on a Minneapolis Honeywell "Brown" Recorder. The 
recorded curve is the experimental cooling curve from which the maximum 
slope was determined. 

The mass rate of flow of air was measured by an ASME orifice meter 
with pressure taps one diameter upstream and one - half diameter down- 
stream, using a range of orifice plates. 

Flow friction performance of the matrix was obtained by measuring the 
isothermal pressure drop across the matrix and the mass rate of flow of 
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air through the matrix. Pressure taps upstream of the matrix and upstream 
of the orifice meter were used to record static pressure entering the 
system and at the orifice section respectively. 
Data reduction relations are discussed in APPENDIX A. A sample 
calculation for one test run is presented in APPENDIX D. 
The following data was recorded for each test run; 
(a) Orifice diameter (d.) 
(b) Entering static pressure (P.) 
(c) Pressure drop across the matrix GA5r=) 
(d) Static pressure upstream of orifice (P) 
(e) Pressure drop across the orifice (AP) 
(f) Orifice temperature (t.) 
(g) Atmospheric pressure (P tm 
(h) Time - temperature cooling curve 


(1) Recorder chart speed 
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4. DESCRIPTION OF SURFACES 

Surfaces examined in this experiment fell into one of two major classi- 
fications - perforated nickel plate matrices or offset rectangular fin of 
stainless steel. In the former case, two of the specimens examined - 
the 20 ° skew matrix and a modified parallel plate matrix, both of 160/40TV 
perforated nickel - had previously been examined by Bannon [1] . These parti- 
cular surfaces were selected to provide an additional check on the system 
reliability. 

The 20 ° skew matrix was formed by corrugated sheets, which were skewed 
relative to one another to prevent meshing of adjacent sheets at an angle 
of 20°. This matrix configuration unlike both the modified parallel 
plate and the offset rectangular fin matrices does not require splitter 
plates. | 

The "modified parallel plate" configuration consists of formed sheets 
separated by splitter plates to prevent nesting. The assembled matrix is 
a maze of rectangular passages with a value of < , the aspect ratio of a 
rectangular flow passage, = 7. The following sketch depicts the ideally 


formed "modified parallel plate" matrix: 





A modified parallel plate matrix of each type of perforated nickel was 
constructed and tested to examine the effects of varying shape, quantity and 
size of perforations on this configuration. 

The offset rectangular fin surfaces were provided by the courtesy of 
the AiRsearch Manufacturing Company, Los Angeles. The method of identifi- 


cation employed with this configuration is best illustrated by the follow- 


ing example: [2] 


16R ae ~153 °143 (0) = 004 oo SS 


OO @ ® 


Number of fins per inch 

Type of fin flow cross-section (R = rectangular) 

Fin height, inches 

Uninterrupted fin length in the flow direction, inches 
Type of surface ( 0 = offset) 


Fin metal thickness, § , inches 


Q8OOG8 O00 


Fin material (SS = stainless steel, Al - Aluminum) 

Two matrices of this configuration were investigated: - the 16R - 
-153 - .143 (0) - .004 - SS and the 20R - .100 - .125 (0) - .004 - SS. 
The offset - rectangular fin matrix is comprised of offset fins, separated 
by splitter plates (0.004" thick) to prevent nesting of the fins which 
would alter the flow characteristics. The isometric sketch below best 


identifies this plate fin arrangement: 


-_ 


SES 
a Ll 
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Dye PRESENTATION OF RESULTS 

Heat transfer and flow friction data for each matrix investigated are 
exhibited in both tabular and graphical form, TABLES I through X and 
Figures 9 through 18, using the Colburn j modulus, Fanning friction factor, 
and Reynolds number. In evaluating the Reynolds number, the hydraulic 
diameter, which is defined as four times the hydraulic radius, was used. 
The effects of entrance, exit, and flow acceleration have been considered 
in the evaluation of f. (See Appendix A). 

Geometric and physical properties of the surfaces are given in Figures 
3 through 8. Figures 19 through 23 show heat transfer and friction compari- 
sons between various matrices investigated. In Figures 24 and 25, "figure 
of merit'' type curves are presented. The j/f versus Reynolds number thus 
plotted gives an indication of the required matrix flow frontal area for 
a given pressure drop. 

Lastly, in Figures 26 and 27, heat transfer power versus flow friction 
power curves on a unit area basis, evaluated for fluid properties at stand- 
ard conditions of dry air at 500 °F and one atmosphere, are represented. 


(See Appendix A). 
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6. EXPERIMENTAL UNCERTAINTIES 

Those items of possible error arising from lack of exact duplication 
of the idealizations, boundary and initial conditions are discussed in 
the "Experimental Technique" portion of this report and no further attempt 
was made to evaluate this type of error. 

The experimental errors which do not fall into the above area are 
three fold and are discussed below. These three categories are: (1) 
uncertainty of physical constants, (2) inaccuracy in the geometric meas- 
urements and (3) instrumentation inaccuracies. For the determination of 
the uncertainties reported, the method of Kline and McClintock [12] was 
used. 

(1) Values for the physical constants required were obtained 
from (4 ] ; [5] and [6 | - The uncertainties in these values, as best 


as can be determined, are listed below: 


co ot 0.5% 
Cc + 0,27 
Pp — 
Np. + 2.0% 
k + 0.5% 
oe 
1.07% 


= 
I+ 


It is noteworthy that the accuracy for the value of C. depends upon 
the metal used. Where pure metals are used the accuracy quoted is feasible, 
however, this is not the case with alloys. Thus the tolerance of c. can be 
a major cause for inaccuracy [13] : 

(2) Due to inconsistencies in construction of the matrices, 
errors in lineal dimensions are considered to be + 0.5%. Inasmuch as the 
weight of the matrix can be determined as accurately as desired, the error 
in W, is considered negligible. Based on this, the uncertainties in geo- 


metric measurements are considered to be as follows: 
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A, Agys “ A_, Ay, + 1.0% 
L + O56 
" < 0.1%, negligible 
(3) Instrumentation inaccuracies were essentially those 
obtained in pressure measurement. With the exception of the temperature 
measurement at the orifice, all temperatures recorded were temperature 
differences which were recorded in inches (See Appendix A). The tempera- 
ture measurement at the orifice was obtained by a copper-constantan shield- 
ed thermocouple which was read by a Leeds and Northrup portable potentio- 
meter. Assuming adequate manufacturer's calibration of the thermocouple 
wire, the estimated possible error is + one-half of the smallest division 
on the potentiometer, or + 0.0025 millivolts which is approximately + 0.1 
°F. 
Since the range of pressures varied, requiring several different mano- 
meters, the maximum uncertainty encountered was used in the uncertainty 


analysis. These values are listed below: 


Po 1 1.25% 
LP, i leon 
OP nF 1.70% 
mM 
te 
atm a 0.0005" Hg (negligible) 


One final parameter for which uncertainty must be determined before 
an analysis of the experimental results can be made is Ney Inasmuch as 
Colburn j factor is not a linear function of Reynolds number, uncertainty 


for high and low N, values was considered. Using Figure 30 and the un- 





R 
certainty in Maximum Slope as 2.0% these values were determinable. 
dNau 
roe Kia AO oe - = + 7.2% 
Tu Tu 
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d 

For N, = 25.0, A= 0.05 Ntu =+ 10.0% 
fom > 
Tu 


Now, with the aid of reference [12] , the uncertainties in j, m, Ne? and f 


were evaluated based on 20: 1 odds. 


2 ifs 
We = (aR) + (BS) sooo (Bw | 


where w_, Wo» eee ee uncertainty interval based on certain odds, 


and We is the interval for the result. 
Considering the evaluation of the Colburn j factor: 


A 2/3 
is WN => N 
J Tu A Pr 


then: 


oe ip 
3 | Wy d) eva) + (24 Wn . 
ou «)* (3 xv) (Sie m 


= 
i 


7 Ae Or Lee 
(As Nov Wa) + ( Mey Wn) r (Ss Ner ns) + 
A e 


} 
Zs N ee ‘ a 
( = Nru fe Pr Wngt) 


Dividing by j yields: 


2 gh >) ~~ vm 
Wi = wie | ne eS fe ( whe at (2 | 
J Nru A Ac 8 Ner 2 
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For the N —seeeeiry OF nN ae0 cas: 
Tu 


A 
: 2 2 4 ° 
= = [ (ore)+ Cor} +601) . (F) ] aaa J — 
For the N Se 25008 a = 0.05 case: 
Tu 


= . Ke 1) + (01) 4 (.o1) + cay 1018 


By similar analysis the uncertainties at 20:1 odds were determined 


ve 
Ute 
=a 


wi = tow 


for the following: 


m + 1.0% 

+ e2 
No + 2.3%, 
f + 4.3% 
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7. DISCUSSION OF RESULTS 





In order to verify the reliability of the data obtained by the test 
apparatus used, several surfaces on which data was currently available 
were tested. Two of these matrices, the 160/40TV perforated nickel 
"parallel plate" and 20° skew, were evaluated by Bannon [1] . These two 
matrices were originally evaluated by the same transient test technique 
on basically the same test apparatus; therefore, satisfactory comparison 
of the data provided the necessary confidence in the test rig. Figures 
9 and 16, as well as TABLES I and VIII, present the results. 

The two offset rectangular fin matrices tested were essentially tested 
by Briggs in [2] with the exception of slightly thicker splitter material 
and using aluminum in lieu of the stainless steel. The data in [2] was ob- 
tained by a steady state "steam-to-air" test technique. Although this 
technique permitted investigation at#much higher values of Reynolds number, 
a sufficient overlap in the Reynold's range existed so that a comparison 
was possible. A good comparison for the 20R -,100 - .125 - .004 - SS 
surface for Colburn j factor exists; however, the friction factor re- 
ported in (2 ] differs noticeably. Results are shown in Figure 17 and 
TABLE IX. For the 16R - .153 - .143 - .004 -SS surface, the values for 
Colburn j factor for Reynolds number < 1000 are in excellent agreement 
with those presented in [2] - In the transition range, the j values re- 
ported herein are noticeably Pye however, the values for friction 
factor agree nicely eHEOGeneer the entire duplicated Reynolds number range. 
It is important to note here that by the "maximum slope" method of analysis 
an inherently large error is possible in the low New region (1.5 < es 
3.5) [7] . The size of this error in N relative to an error in maximum 


u 


slope can best be appreciated from Figure 30. Consider, for example, a 
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2% error in maximum slope with an Neu value of 2.3 and AX O. The 
percent error in See from Figure 30, is of the order of 20%. For values 
of Nou between 2.0 and 2.3, the percent error has not been evaluated. The 
156 Nru < 3.5 region coincides with the Reynolds range for which the 
reported j values differ from [2] » and may account for the difference 
noted. The results for this matrix appear in Figure 18 and TABLE X. 

In spite of the areas of disagreement, the correlation was such that 
the data for the j in the Reynolds range < 1000 and the friction factor 
data for these matrices are considered reliable. It is noteworthy that 
the surface tested in [2] is specified as 15.76R - .153 - .143 -.004 - Al, 
so that the matrix evaluated herein may only be an approximate duplicate. 
A comparison of the relative performance for these matrices is presented 
in Figure 23. Additionally, Figure 25 illustrates j/f vs. Np» a flow 
area "goodness" factor. 

The perforated nickel matrices investigated covered a wide range of 
percent open area, 12 to 50%. Both slotted and round perforated nickel 
surfaces were evaluated for the "modified parallel plate" geometry. The 
perforated nickel geometric properties are given in APPENDIX C for the 
various types investigated. In order to make a reasonable comparison of 
the performance among these matrices, the same frontal area and flow length 
were specified. Results are presented in Figures 9 and 10 and listed in 
TABLES I and II for the slotted perforated nickel matrices; those for the 
round perforated matrices are given in Figures 11, 12, and 13 and tabu- 
lated in TABLES III, IV, and V. 

In order to establish a meaningful fiducial point, a solid nickel 


matrix of the identical material (obtained from the same manufacturer) 


with the same frontal area and flow length as the perforated matrices was 
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fabricated. The performance of this matrix is illustrated in Figure 15 
with tabular values given in TABLE VII. 

The best performing slotted perforated nickel matrix was the 160/40Q. 
The 50G round perforated nickel (50% open area) was the best surface in- 
vestigated among the "parallel plate" nickel matrices. 

As another approach to evaluating the Colburn j factor, j' was con- 
sidered. The parameter oe based on the plane surface area, A*, as if 
no perforations were present; whereas, j is based on the heat transfer sur- 
face area, A, which includes the solidity correction, s. The value of j' 


is readily attainable from the relationship 


j' =jxs , where s = A/A*. (See Appendix D). 
The determination of the solidity is fully described in APPENDIX C. Figure 
13 for the 50G perforated nickel matrix includes the values for j' versus 
Reynold's number in addition to the f and j values. 

A comparison of the performance of the two slotted perforated matrices 
together with the solid nickel matrix is shown in Figure 19. A similar 
comparison for the round perforated nickel matrices is given in Figure 20. 

It is readily apparent that the perforated material yields a higher 
heat transfer characteristic. Inasmuch as the friction factor also shows 
an increase, a more appropriate comparison can be made from the flow area 
) versus 


"soodness" factor curve in Figure 24. The heat transfer power (h 


STD 


flow friction power (E.,,.) per unit area given in Figures 26 and 27, for the 


STD 
slotted and round perforated nickel matrices respectively, is another ''good- 
ness" factor evaluation. Here the information for each of the matrices is 
presented based on standard conditions for fluid properties and a common 


hydraulic diameter. (See APPENDIX A). 
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So that an appreciation for the perforated material might be realized 
in direct transfer type of heat exchanger application in addition to its 
use in matrix type for rotary regenerators, a "parallel plate" matrix with 
the 160/40TV perforated nickel fins and solid nickel splitter plates was 
fabricated. The performance of this matrix is given in graphical and 
tabular form in Figure 14 and TABLE VI. A comparison with the solid 
nickel matrix is made in Figure 21. Figure 24 contains the j/f presenta- 
tion and Figure 26 shows the er vs. E orp for this matrix. Here again is 
a noticeable increase in heat transfer performance, but with a negligible 
increase in friction factor. 

For the 160/40TV perforated nickel material, a 20° skew matrix was 
compared in performance against the "modified parallel plate" to see how 
these two configurations compare with one another. The separate perform- 
ance for each has been mentioned previously; however, the direct comparison 
is shown in Figure 22. As with the other nickel surfaces, the flow area 
"goodness" factor is represented in Figure 24. 

In the results of this investigation are several items which are general 
in nature, but nonetheless extremely important. In the laminar Reynolds 
number range, friction factor is inversely proportional to Reynolds nunber. 
By Reynolds analogy, it was anticipated that the Colburn j factor should 
likewise be inversely proportional to Reynolds number in this region. 
However, the experimentally obtained values of j in the low laminar Reynolds 
range utilizing the equivalent conduction length (defined in APPENDIX C) 
together with TABLE XI or Figure 29 yielded lower values for j. For this 
reason, a re-evaluation of the conduction parameter, a and maximum slope 
data points for Figure 29 was performed. The equivalent conduction para- 
meter AL was computed and the Nou value was determined by utilization of 


. This new 


the extrapolated j vs. Ne curve, since j is proportional to Neu 
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determined value of Neu was considered to be the actual Ba, for the parti- 
cular "maximum slope" and conduction parameter. Values determined in this 
fashion are asterisked in TABLES I through VI and VIII. The tabulated 
conduction parameter values in these tables are the equivalent conduction 
parameter (See APPENDIX D). 

The second item is that of an equivalent conduction length which was 
quoted earlier. In the perforated fin material the slots or round holes, 
as the case may be, provide for a winding heat flow path. It was appreci- 
ated that the matrix flow length, L, was not the true conduction path but 
somewhat less than actual. A mean conduction path length evaluation is 
presented in APPENDIX C. The conduction parameter, A. , based on the flow 
length, L, was evaluated, and by simply multiplying by the ratio L/L (also 
given in APPENDIX C), the value of hi was determined. 

A third item is the inherent error in Ne in the range 1.5< N i= ce) 


T 


which was mentioned earlier. Howard [7] published an error curve for N_ , 
Tu 


including slope vs. Nes see Figure 30. 

Fourthly, the entrance and exit length effects for a matrix with a 
small hydraulic radius are small. In the evaluation for Fanning friction 
factor (APPENDIX A), the effects of exit and entrance length were consider- 
ed using values of K. and K, from Figures 37, 38, or 39, as appropriate. 
For the perforated surfaces, the K values were taken corresponding to 
laminar flow conditions; however, in the offset rectangular fin geometry 
the values were taken from those Figures for BE = OO in accordance with 
reference [11] . The offset fins never have velocity profiles that are 
fully established, that being a purpose in the design. The same calcula- 


tions were repeated on the offset rectangular fins using the K values for 


laminar flow. The friction curves plotted for the offset rectangular fins 
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were coincident with those obtained earlier; however, some change was noted 
in tabulated values. 

A final item worth mentioning is again concerned with the low Reynolds 
range. A slight alteration in the conduction parameter value has a marked 
effect on the a value in this range, which is apparent from Figure 29. 
For this reason, when utilizing a nickel material, or any other material 
with a wide range of specified values for thermal conductivity, it is 
mandatory that the appropriate value of Ke be obtained for calculating 


purposes. The value of A is directly proportional to this k, value, 


i.e., 
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8. CONCLUSIONS 

Reliable data is obtainable from the experimental apparatus based on 
a comparison with those surfaces previously evaluated elsewhere. 

A very definite improvement in performance is achieved with the utiliza- 
tion of a perforated material. Kays, in reference [10] » postulated the : 
improved heat transfer characteristic without a "large amount of form drag 
so characteristic of high performance surfaces". The perforations apparent- 
ly disturb the heat transfer sub boundary layer so that tubulence increases 
and the resulting high convection conductance is realized. Inasmuch as the 
perforations are fine, the downstream edge is considered to be in the wake 
of the separated region. The perforated nickel matrices investigated here- 
in confirm Kays hypothesis. Unfortunately, inasmuch as the perforated sur- 
face investigated in [10] is greatly different in geometric and physical 
properties, no meaningful comparison with those tested in this investiga- 
tion could be made. 

As stated previously, the experimentally obtained Colburn j values in 
the low laminar Reynolds range using TABLE XI or Figure 29 were lower than 
j values obtained by extrapoldeinc the j versus N, curve. Inasmuch as fric- 
tion factor and Colburn j factor are both inversely proportional to Reynolds 
number in the laminar region owing to Reynolds analogy, the values for 
Colburn j factor extrapolated from the experimental curves are considered 
more reliable in the low laminar range. 

Effects of longitudinal conduction have been considered in the analysis 
of data. For the perforated material matrices, corrections for the actual 
conduction path length and for the area removed by the perforations have 
been included. 

The best performance for a perforated "parallel plate" matrix was at- 
tained by the 50G type material. For the offset rectangular fin surfaces, 
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ix was superior. 
- .125 - .004 - SS matr 
the 20R - .100 
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9 RECOMMENDATIONS FOR FURTHER STUDY 

1. Due to the inherent error in the high Reynolds number range where 
1.5 Nu = 3.5, a more reliable means for determining Neu values should 
be investigated. Perhaps pursuit of Hausen's cyclic operation [9] , a 
steady state periodic-flow technique, is appropriate in that this method 
measures effectiveness from which 2+ can be determined subsequently. 
Since the results obtained using Schumann's analysis, even with longitudin- 
al conduction effects considered, at the low Reynold's range appears in- 
adequate; the use of another testing method may prove fruitful. 

2. Inasmuch as the best performing surface for perforated nickel 
was the 50G type (50% open area), other types, such as the 40G which has 


a slightly larger hole diameter, with same plate thickness and percent 


open area as the 50G may prove worthy of test. 


26 





10, BIBLIOGRAPHY 


{1} Bannon, John M. 


"An Experimental Determination of Heat Transfer 
and Flow Friction Characteristics of Perforated 
Material for Compact Heat Exchanger Surfaces". 
Unpublished Master's thesis, U. S. Naval Post- 
graduate School, Monterey, 1964. 


[2] Briggs, D. C. and London, A. L. "The Heat Transfer and Flow Friction 


[3] Coppage, J. E. 


[4] Eldridge, E. A. and 


Characteristics of Five Offset Rectangular and 
Six Plain Triangular Plate-Fin Heat Transfer 
Surfaces'', International Developments in Heat 
Transfer, Part I, International Heat Transfer 
Conference, Boulder, Colorado, 1961. 


“Heat Transfer and Flow Friction Characteristics 
of Porous Media", TR. No. 16, Department of 
Mechanical Engineering, Stanford University, 
Stanford, Calif., December 1, 1952. 


Deem, H. W. Report on Physical Properties of 


Metals and Alloys from Cryogenic to Elevated 
Temperatures. Philadelphia: American Society 


for Testing Materials, 1961. 


[5] Goldsmith, A., et.al. Handbook of Thermophysical Properties of 


Solid Materials. Revised Edition, Vol. I, 
Elements, New York: The MacMillan Company, 
1961. 


[6] Hilsenrath, J., et. al. Tables of Thermodynamic and Transport 


[7] Howard, C. P. 


(8) Howard C. P. 


[9] Jakob, M. 


[10] Kays, W. M. 


Properties of Air, Argon, Carbon Monoxide, Hydrogen 
Nitrogen, Oxygen, and Steam. New York: 


Permagon Press, 1960. 


"The Single Blow Problem Including the Effects 
of Longitudinal Conduction'', ASME paper, No. 
64 - GTP - ll. 


Bureau of Ships, Code 645, U. S. Navy, private 
communication, February, 1965. 


Heat Transfer, Vol. II, New York: John Wiley 
and Sons, Inc., 1957. 


"The Heat Transfer and Flow Friction Characteristics 
of a Wavy Fin, a Strip Fin, and a Perforated Fin 
Heat Transfer Surface", TR. No. 39, Department of 
Mechanical Engineering, Stanford University, 
Stanford, Calif., October 31, 1958. 


27 





[11] 


[12] 


[13] 


[14] 


[15] 


[16] 


[17] 


[18] 


[19] 


[20] 


Kays, W. M. and London, A. L. Compact Heat Exchangers. Second 
Edition. New York: McGraw-Hill Book Company, 
Inc., 1964. 


Kline, S. J. and McClintock, F. A. "Describing Uncertainties in 
Single-Sample Experiments", Mechanical 


Engineering, January, 1953, pp 3-8. 


Locke, G. L. "Heat Transfer and Flow Friction Characteristics 
of Porous Solids'', TR. No. 10, Department of 
Mechanical Engineering, Stanford University, 
Stanford, Calif., June 1, 1950. 


London, A. L. "New Developments in Compact Exchangers - Design 
Theory, Surfaces, and Applications", lecture 
presented at 6th National ASME - AIChE Heat 
Transfer Conference, Boston, Massachusetts, 
August, 1963. 


Murdock, J. W. "Tables for the Interpolation and Extrapolation 
of ASME Coefficients for Square-Edged Concentric 
Orifices", ASME paper, No. 64 - WA/EM - 6. 


Power Test Code Supplements, Instruments and Apparatus, PTC - 19. 
5.4 Flow Measurement, Chapter 4. New York: 


American Society of Mechanical Engineers, 1959. 


Schumann, T. E. W. "Heat Transfer: A Liquid Flowing Through a 
Porous Prism", Journal of the Franklin 
Institute, Vol. 208, 1929, pp. 405 - 416. 


Smith, R. H. and Wang, C. T. “Contracting Cones Giving Uniform 
Throat Speeds", Journal of the Aeronautical 
Sciences, October, 1944, pp. 356 - 360. 


Solar Aircraft Company, Engineering Report ER 1221, Regenerator 
Core Module Tests - T - 600 Engine, 18 April 


1962. Rev. A: 25 October, 1962, AD 290234. 


Theoclitus, G. and Eckrich, T. L. "An Experimental Technique 
for Determining the Effectiveness of Extended 
Surfaces'', presented at 7th National Heat 
Transfer Conference, Cleveland, Ohio, August, 
1964. 


28 





snjzerteddy ysay jo wersetg ot3ewWeYUdS 





T € t) i 
iL oy L a9 / 
_—<.<— ° °° © q 
/ 
| / 
| 
: | bas: 
; | 
] 
_ Ff 
z 4 
A9UdAYSTeAAS ie d P 
use220¢5 P 
eo \ uoT30e6 
‘ quaweanseay Mold 


sde], 
sanssaig XTAqeW 





*[ sangty 


DATEA 
[or}u0D 
MOTTA 


\ 
| 


BOTFTAO 


a as 


PATPA RETUI 
atossoadmoy 








WO © 





LEGEND 


Honeywell Temperature Recorder 
W20H Variac Autotransformer 
Orifice Metering Section 
Precision Potentiometer 
Heater Selector Panel 

Orifice Thermocouple 
hianometer Board 

Thermocouple Selector Panel 
Test Section 


See et ne 











afin s 

a¥a AT YI 2 
~. * 
wpe we ok 

a a ve 
snermna 
+f, 

dba ch. | 

NT, 

LOU 

“oF 

¥ ‘ 

» ye re 
rF 4 
aBd Nee r + A 
VW wr x 
abe _ - 


DI 


 o 


hw 





f f 


LOU / é,' a¥ 1OU/ UU 


Perf.Nict Perf Nickel 


e etd / 


‘alee 
x 
e e 
. e ad & 
> 
e “ e 
bee 
e . 6 4 
, s 





€88°0 
8°996 
820200 °0 
9°TE 

; 885 TT0°0O 
Z9T90°0 
€7Z00°0 
€5690°0 
€€0z° TT 
9T00°0 
L°s¢ 
S90T°0 


T9X°FN “Sted 
D 0S 


SOoTszeW yOd®Td TeTresedy ( 90S pue *d Set *W SZT ) TeXx°TN 


peyeIojiseg punoy jo setysedoig [eotshyg pue ofFIyeUOen °4Y amnsTy 


7$8°0 

L° Lest 
T96100°O 
eT 
88STTIO°O 
6£65 0°O 
07900°0 
€5690°0 
%c0L° LT 
0z00°0 
“LS 
$90T%0 


Tex°TN * sated 
d Scr 


€88°0 

0° €O7T 
8Z0200°0O 
wt 
88STIO°O 
ZT190°O 
G€700°0 
€5690°0 
SLSe°9T 
9T00°O 
L°8€ 
S90T°O 


TSXOTN *jsed 
W S2T 


(d) A4}s010g 


¢33/>35 (9) sseuzoedmog) 


43. (Hq) seqowetq oF Tnespsy 


cw /at ("38) fa;sueg xPrIeR 


cu (TA) SUNTOA XFIIEF 


oi3 (OV) Bory MoTY Seay 


o13 (ly) very uotqonpuog TeI0] 


zw (“y) Beary Teqoag 


o33 (¥) Peary Jessuery JOH TBO] 


soayouytT *‘ssouyoTuL [Teyse.ey 


Jo WU /naq (84) AaFATZONpUCD TemIeUY 


Jo AT/n4g (59) qeaH oTFFOEdg 


Tepseqzey XFS 


~~ 


32 








7N\ 








at, ~ + } 


Matrix 


Spe cific 


Thermal BPOmMeaueOULY Luy 9) 
(25 + 1 " - 
Material KNESss, 


Jeera) 6 i iv i 


> 


Free Flow ea (A,) 
Matrix Volume ({ 
Matrix Density ( ¢ 
Hydraulic Dj 
Compact ( 


Porosity (p) 


fi aD | ° 


¥ 








| 


= | i ba | 

i 2 ames 7 ee 3-4 
THE L.S.STARRETT Co, 3 
| 


f= 





20 Degree Skew 


160/40 TV 
Matrix Material Ferf Nickel 
Specific Heat (c,) Btuylip 0.1065 
Thermal Conductivity (kg) Btu/ hr ft F 38.7 
Material Thickness, inches 0.0022 
Total Heat Transfer Area (A) ft* 15 del 
Frontal Area (Ap,) ft” 0.0734 
Total Conduction Area (A, ) ft< 0.00349 
Free Flow Area (A,) ft 0.06412 
Matrix Volume (V,) ft? 0.01211 
Matrix Density (@,) lb/ft? 505m 
Hydraulic Diameter (Du) ft 0.00252 
Compactness (®@) ft2/rt- 1272.8 
Porosity (p) 0.873 


Figure 7. Geometric and Physical Preperties of a Slotted 
Perforated Nickel (160/40 TV) 20° Skew Matrix 
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Surface Heat Transfer and Friction Data 


Figure ll. 
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Surface Heat Transfer and Friction Data 


Figure 23 


"Parallel Plate", 125 P Perforated Nickel 
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Figure 15. Surface Heat Transfer and Friction Data 
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Figure 16, Surface Heat Transfer and Friction Data — 
20° Skew, 160/40 TV Perforated Nickel 
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Figure 31. Matrix Holder 
(Upstream View) 

















Figure 33. Velocity Profile and Temperature Distribution 
Measurement Device = 
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Figure 34. Nichrome wire Heater 
(Upstre+m View) 
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Close-Up of Apparatus Test 
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Figure 37. Entrance and Exit Pressure Loss Coefficients 
for Triangular Cross Section Matrix with 
Abrupt-Contraction Entrance and aE 


Expansion Exit 
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Figure 38. Entrance and Exit Pressure Loss Coefficients 
for Parallel Plate Cross Section Matrix with 
Abrupt-Contraction Entrance and Abrupt- 
Expansion Exit 
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Figure:39, Entrance and Exit Pressure Loss Coefficients 
for Multiple - Rectangular Cross Section Matrix 
with Abrupt-Contraction Entrance and Abrupt- 
Expansion Exit 
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APPENDIX A 
DATA REDUCTION RELATIONSHIPS 

Summarized herein are those data reduction relationships which are 
of prime importance in obtaining the heat transfer and flow friction 
characteristics as presented from the raw data taken. Before the data 
on the matrices tested can be reduced, or a comparison of the different 
matrices can be made, it is essential to obtain accurate geometrical dimen- 

| 


Sions. The geometrical factors of interest are the porosity, p, the hy- 


draulic diameter, D,., and the area compactness, GC . Determination of 


H’ 


any of these fixes the third since it may be shown that 


= 4 
C _ [14] (A-1) 


The definitions of these terms are: 


Free Flow Area he 
(1) p = Frontal Area = Aer (A-2) 


4 x Free Flow Area 








(2) Dy = Heat Transfer Surface Area = 4 Ac L (A-3) 
A 
Heat Transfer Surface Area = A 

(3) € = Frontal Area x Flow Length Ag L (A-4) 


MAXIMUM SLOPE 


As presented earlier, the maximum slope of the generalized cooling 


curve is a unique function of Neu and XL. 


(A-5) 


YD (Nru A) 








7T and N,, have both been previously defined, but for convenience 


the relationships are restated: 





ee) hA S, and N+ = h A 


Wy Cg MCp 





Therefore: 


TAN Be (ri Ce/Wece | © 


and d(7/Nru) = mep/We ce) de is (c/é, ) Ales (A-6) 
Furthermore: 


tp.7 ti 1G Te 
2) I ASFap es + | 7) 


and 





- 1 __ diag -t 
eee Et Alaa [te a | (A-7) 


Combining (A-6) and (A-7) 


teg- ti = 3 
ae ta ~ \> a a ee tey 
dT / Nu | max 


(A-8) 


From Figure 28: 


Se d{te2- ter | ay 
dO ae . 
(2) = dO » sec 


a 
Ri) 


eer 
Chart speed 





(3) y=d(te2-tyy ), inches 


(4) M= (te, - ti), inches 


Combining these determined values into (A-8) yields: 








toy -t; = 
Cy : i ,» . Chart speed (A-9) 
, dl y/ Neu | G M . 
| imax 
With this value of maximum slope and A_ = Ka Ps , enter TABLE 
p 


XI or Figure 29 to obtain a corresponding SF value. 
u 


MASS RATE OF FLOW: 
The mass rate of flow is calculated from the method given in reference 
[16] with necessary rearrangements in order to reduce unnecessary hand 


calculations and to gain the full benefits of Murdock (15 | : 


359Kd, * F. ¥ \ AR ¥ (A-10) 


m = 
where: K = flow coefficient, including velocity of 


cS 
approach = |—>—____ with C = orifice 
Vi- G4 


coefficient of discharge (obtained from [15] ) 
and 6 = ratio of orifice diameter to internal 
pipe diameter, d /d 

d 


0 


F 
a 


orifice diameter, inches 


thermal expansion factor 
Y = expansion factor 


P 
¥ = specific weight of fluid flowing = RT 


2 
with P in #/ft 
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R = 53.35 ft - #/lb °R for air 
T = t(°F) + 459.7 = degrees Rankine (°R) 
Exes = Pressure drop across the orifice, inches H,0 
substituting — for ¥ » since the magnitude of the local acceleration of 


gravity, g, is taken as being equal to the magnitude of standard accelera- 


C for K 
{ ~ Q* 


tion of gravity, g 3 and 


gives: 





.. 359C d’ FE Y\AR EP 
mM — 'e) a 
vil ieF ial (A-11) 


Inserting the numerical value for R and converting P to units of 


Pita’ yields: 





= 589. BL CU FLY 


me Ye 


From Figure 38, reference 16] : F~1.00 


A EKpsie) (A-12) 





From Figure 40b, reference [16] : Ye 0.9985 ~ 1.0 


The pressure, P, in (A-12), is hereafter referred to as P, and is 


determined by the following relationship: 


Py Fi Pen 


Where P atm is the atmospheric pressure in inches of mercury, 


~ P/ 13.6) 0.4912, #/in? (A-13) 


ae is the static pressure upstream of the orifice in 


inches of H,0 
13.6 inches of H0/ inches of mercury 


0.4912 #/in?/ inches of mercury 
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The temperature, T, in (A-12) is the air temperature at the orifice 


in degrees Rankine (T,): 


Le cae 459.7 (A-14) 


Replacing re and Y in (A-12) with their numerical values, the equa- 


tion reduces to the working form: 





m = (A-15) 
Per 
REYNOLDS NUMBER 
Reynolds number is defined as: 
Z Si 
Neg oat (A-16) 


ye 


where G is the mass flow velocity based on the free flow area, A.- 





ae 
G Ae Pp Aer (A-17) 


Substituting (A-17) into (A-16): 


N = a Dy 
i Per Li 


(A-18) 
me : 1 
But from (A-1) D/P = G , therefore: Ne = ae (4) (A-19) 
PAec\ @ 
Equation (A-19) shows that for a given matrix and mass flow rate the 
Reynolds number is inversely proportional to the compactness. Substitut- 


ing (A-4) into (A-19) yields the working equation: 


= 4Lm 


7 ae (A-20) 


Nr 





FANNING FRICTION FACTOR: 


The following sketch together with equation (A-21) describes the flow 


system under consideration: [11| 





2 V Pas Vin 
bP = aor van (K Lar = VA ie ea oe 1) So - A. ail 
entrance flow core 
effect acceleration Pprepion 
vO. 
- (1 ~ p® - Ke) VY} (A-21) 
exit 
effect 


It is important to note that Ves v1 and VY, ~ Me in the above sketch 
as the pressure changes (for gas flow heat exchanger application) from 
section 1 to a and b to 2, respectively, are very small relative to the 
total pressure. Inasmuch as the testing is performed with air at moder- 


ate temperatures and pressures, the perfect gas law is applicable (Peony. 


Substituting v = Ne and recalling that oT = Ac or Cy - Ac 
. 4L A C 


from equation (A-3), equation (A-21) solved for f for the isothermal 
case becomes: 


AP |. Ea & + -) 
f 2 cin oS ee ; Fy Pe 


con’t on 
next page 
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(con’t.) 


( + Pa ( 2 2 2) (a ‘ ) | ry 
2 Po en a (A-22) 


where Qe m = Wald with the upstream values subscripted by 1 
and the downstream subscripted by 2. K, is the entrance coefficient and 
K, the exit coefficient, and both are dependent on porosity, shape of 
the flow cross section and the matrix Reynolds number. Values a K. and 
K, are obtained from Figures 37, 38, or 39, which were reprinted from 
reference [12] : 

If we consider an order of magnitude approximation, the first term in 
(A-22) is by far the greatest contributor to friction factor for small 
pressure differentials. 

The approximation, therefore, that P+P, _ pee b 


reduces (A-22) to the following 


1s =[2ec On( SE) = (Ko + Ke )+ BE (1 + pa) ae (A-23) 
L 


Reference 114 | points out that the relationships given here can be 
recombined in such a manner so as to show that f is proportional to the 


porosity cubed and inversely proportional to the compactness cf OC p?/ @ Dee 


COLBURN J FACTOR: 


Colburn j factor is defined as: 


2 
= SL Ney . (A-24) 
GCp 


2/ 
j= Noy Nor : 
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Substituting (A-17) for G and multiplying by A/A yields 





=a HA Ae 2/3 
j 7 Np, (A-25) 
but LA 
Cp = 2g 
therefore: 
; 2 
j= Nero = Ner ‘3 (A- 26) 


Combining equations (A-1) and (A-3): 


thus: 


P (A-27) 


It is apparent from equation (A-27) that the Colburn j factor is 


directly proportional to porosity and inversely proportional to compact- 


ness Cj of P/Q) 


HEAT TRANSFER POWER AND FLOW FRICTION POWER: 


An evaluation of the heat transfer power versus flow friction power 
is of interest in that it is a measure of relative performance. 


The heat transfer power per unit area per degree temperature differ- 


ence is Ju | : 


h= C (A-28) 
bal 7c Se? nN 
N 4/3 Du J R 


Pr 
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Equation (A-28) with the properties evaluated at standard conditions 


of dry air at 500 F and one atmosphere becomes: 


Deny = 0.02195 (5. )(Nr j) Beu/br ft? °F (A-29) 
where > 0.2477 Btu/1b °F 
P= 0.0678 1b/hr ft 
© = 0.0413 1b/fe? 
N, = 0.671 


The flow friction power per unit area is [11] : 





3 
1 ee (+) c Ne (A-30) 
2 &c ia Du : 


When equation (A-30) is evaluated at standard conditions oe and e 
values given above), the flow friction power per unit area at standard 


conditions is: 


Oe xe 0s +) (Ne) HP/Et 
0) mn Du 1000 ‘ 


For comparison purposes the surface geometries were reduced to a 
common hydraulic diameter of Dy = 2x 107? ft. Figures 26 and 27 show such 


comparisons for the nickel surfaces. 
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APPENDIX B 
DESCRIPTION OF EQUIPMENT 

Heat transfer data was obtained utilizing the "single blow" technique 
which consists briefly of monitoring the fluid temperature response at 
the test matrix exit while subjecting the inlet to the test matrix to a 
step change in fluid temperature. 

Friction factor was determined from pressure drop data obtained from 
static pressure taps located in the test section - one at inlet, the other 
at the exit. 

Necessary equipment to perform such an experiment falls into one of the 
following categories: 

(1) Fluid source 

(2) Flow metering system 

(3) Temperature measuring system 
(4) Pressure measuring system 
(5) Fluid heater system 


(6) Matrix holder and test section casing 


FLUID SOURCE: 
Air, the working fluid, was provided to the test apparatus by connect- 
ing the rig to the inlet of a 30HP, multistage, Spencer Turbo-Compressor 


rated for 550cfm operating on a 220 volt a.c. power supply. 


FLOW METERING SYSTEM: (See Figure 1) 
The flow metering system consisted of an ASME standard orifice section 
constructed for d and d/2 pressure taps in a 3.08" inside diameter metal tube. 


fie]. A wide flow range was obtained by utilizing thin plate concentric 


- 
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orifices of throat diameter 2.310, 1.971, 1.232, 1.540, 0.462, and 0.308 
inches respectively. 

Control of the pressure drop across the orifice was maintained by a 
gate valve downstream of the apparatus, ahead of the compressor inlet and 
another gate valve on the compressor inlet. The blast gate on the compressor 
discharge was pre-set so that the unit could not be operated beyond the full 
load rating of the turbo. 

This turbo-compressor had many distinct advantages over the one used in 
reference [1] : some of the more pertinent being (1) a 550 cfm capacity; 
(2) essentially constant pressure throughout the range of the machine; 

(3) the load decreases as less air is used so that the problem of stalling 


at low flow rates is non-existent. 


TEMPERATURE MEASURING SYSTEM: (See Figure 1) 

The majority of the temperatures were measured with iron-constantan 
thermocouples. These thermocouples use iron for the positive conductor 
and constantan for the negative conductor. Measuring T is a group of five 
thermocouples in series, bound together and insulated from each other by 
teflon tape and inserted in a 1/8" diameter aluminum tube to serve as a 
radiation shield. A small. aperature was cut in the inserted tube facing 
upstream. T, and Ty, are five wire thermocouple grids, connected in series, 
mounted permanently in the test section casing. ae is a movable five wire 
series connected thermocouple grid located in the matrix holder so that it 
is adjacent to the downstream side of the matrix. The purpose of connecting 
the thermocouples in series was to magnify the emf output so that the instru- 


ment sensitivity to small temperature changes would be enhanced. With T, 


bucked against T,, it is possible to determine the temperature uniformity 


89 





across the matrix. t versus T, indicates the difference between upstream 
and downstream temperatures; thus, it is this combination which is used to 
record the rate of change of upstream and downstream temperatures with time, 
resulting from the step change in fluid temperature upstream. This differ- 
ential is recorded on a Minneapolis-Honeywell "Brown" Recorder and serves 

as the primary data for heat transfer evaluation of the matrix. The re- 
corder has variable chart speeds so that acceptable cooling curves for 
maximum slope determination can be generated. This instrument has a span 
adjustment which permits the span to be varied continuosuly from the narrow- 
est to the widest span desired from 0 to 55 millivolts. For the work perforn- 
ed herein, the recorder was precalibrated for a 0 - 3 millivolt scale. This 
unit also has adjustable supression and damping adjustment features. The 
damping adjustment provides a filter network for removing a.c. strays, and 
provides correct damping so that neither undershoot nor overshoot in the 
response curves is experienced. 

A copper constantan thermocouple was inserted into the duct just a- 
head of the orifice to measure ue orifice temperature. This temperature 
was read in millivolts on a Rubicon Company Portable Precision Potentio- 
meter, converted to degrees Fahrenheit and recorded for each run. 

A copper constantan thermocouple grid was inserted in place of ie 
to determine the temperature distribution across the cross section during 
initial check out of the apparatus. This thermocouple arrangement was 
part of a special holder lined with balsa wood to provide a smooth, 


continuous flow passage, which also contained an impact tube used for 


checking the velocity profiles (See Figure 33 ). 


PRESSURE MEASURING SYSTEM: (See Figure 1) 


Pressure taps are located upstream and downstream of the matrix and on 
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either side of the orifice. Each tap is connected to an appropriate mano- 
meter or draft gage via Imperial Company "poly-flo" tubing. Quick closing 
valves were installed at various positions in the line to permit isolation 
of sections and as a safety feature. 
The following instruments were utilized: 

(1) Meriam Instrument Company, 0 - 50" manometer 

(2) Meriam Instrument Company, -8" to +8" manometer 

(3) Ellison Differential Direct Draft Gages, 0 - 6", 0 - 8". 

(4) Ellison Inclined Draft Gages, 0 - 0.5", 0 - 6". 

(5) E. Vernon Hill and Company Type "C" Micromanometer, 

O-= 1.250. 

Any one of these instruments or a combination can be used to measure the 
differential pressures of the orifice or matrix or the required static 
pressures - upstream of matrix and upstream of the orifice. Cross checking 


of the various instruments assured reliable operation. 


FLUID HEATER SYSTEM: (See Figure 1) 

The heater section is comprised of 28 nichrome wire heaters of 0.0031" 
diameter. The heater system was designed to elevate the air temperature 
20 ° F above ambient for a mass rate of flow of 1000 1b/hour. The nichrome 
material was selected because it has high resistivity, low thermal conduc- 
tivity and specific heat, and has a very small time constant so that it 
permits one to approach the idealized step function. [1] 

The wire heaters are wound two to a bakelite frame, 1/32" between each 
wire, with 50 and 52 wires respectively to each heater on a frame. These 
heaters are connected in parallel, then via a switch to a variable voltage 


bus. The heater frame and its electrical connections are apparent in 
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Figures 34 and 35. 
A schematic wiring diagram for one of the two heater complexes is shown 
below (1] - A total of 14 frames are wired in this fashion, thereby put- 


ting all 28 heaters in parallel, permitting two to be switched on or off 


at a time. 


Bus Bar 







in Heater Selector 
Selector Switch 
nel o——_o— 





On-Off 
Switch 






v 
220 Variac 






The number of heaters in use decreases as flow rate decreases. The 
necessary voltage variations are obtainable by a General Radio Company 
Type W20HM "Variac" Autotransformer, 0 - 280 volt, 8 amp load from a 240 
volt 50 - 60 cycle line. 

The recommended procedure is to decrease the number of heaters insofar 
as practical while maintaining a high voltage. When there are only four 
heaters in use, voltage reduction becomes necessary to obtain the required 
20 ° temperature rise, keeping the number of heaters constant. Meeting 
the current requirement is the dictating factor. It is noteworthy that the 
heaters should never be energized unless sufficient air is flowing, other- 


wise the wires will sag and short circuiting will result. 
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MATRIX HOLDER AND TEST SECTION: (See Figures 31, 32, and 35) 

The matrix holder and test section are made of polyethylene plastic. 
The test section casing holds the T, and T, thermocouple grids permanently 
mounted, as well as the upstream and downstream static pressure taps required 
for determining friction factor. All parts were machined to close tolerance 
Co assure a snug fit for the matrix holder and to guarantee good alignment 
of the flow channel through the heater and test sections. The flow channel 
is 3 - 3/16" x 3 - 3/16" and matrices of flow lengths up to 3" may be 
tested. As mentioned previously under Temperature Measuring System, the T, 
thermocouple grid is located in the matrix holder. 

Correct positioning of the test matrices in the flow channel is en- 
sured by using styrofoam plastic inserts. In addition to the aforemention- 
ed systems, one item remains to be considered - the Inlet Cone and Screen 
Straightener Section. The inlet cone was designed using Smith and Wang 
(18] as a guide. Although this reference pertained to contracting cones for 
circular duct, the resultant inlet cone for this rig proved very satisfactory. 
To assist in guaranteeing a uniform velocity at the throat, a stainless 
steel 60 mesh screen was installed (See Figure 1). Space is available in 
the screen straightening section to add additional screens if they should 


be required for flow rates beyond the range of this investigation. 
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APPENDIX C 
PERFORATED NICKEL GEOMETRIC PROPERTIES 

The material used in the neviorated plate matrices wae pure nickel, 
electro-deposited sheet of integral structure Serre cured by Perforated 
Products, Incorporated. Types 160/40TV, 160/40Q; 50G, 125M and 125P 
were utilized in the matrices - the first two of these having slotted 
Openings and the remaining three having round openings. A brief descrip- 
tion of each of the aforementioned types as specified by the manufacturer 
ils set forth in Table C-l. 

Inasmuch as the perforations in all cases were conical, the use of 
the average slot length or average hole diameter was considered a propos. 
Close examination of each of the different types revealed that this 
average value was somewhat less than the dimensions specified by the manu- 
facturer in each case. In addition to the requirements for determining 
porosity of the plate, a correction for the increased longitudinal conduc- 
tion length due to the perforations as well as a correction to the solid 

conduction cross-sectional area was performed. The following idealized 


patterns for each type used, best illustrate these corrections. 


160/40TV 


oo 


9095" eOL55" 


a 0008" | | 
' 


Mid 





dj 





~0063"! 





Cross Hatched Area, aan = (.0063) (¢(.025) = .0001575 i 

2 
Slotted Area, Any in Ann = (.0008) (.0155) = 0000124 in. 
Plate Porosity, the fraction of open or slotted area, 


A ' 


vy = St 
Sop 
— : A 
Solidity, s, is l-zw =1- sl 
Bra 
.0000124 _ 


a —__—< +) 





0063 nos L; 


= - | oe 


_— 





L , conduction path length, = 0.0125 + 0.0063 = 0.0188" 
2X 
a 0063 


= 0.3351 





Lae 0188 
aL 


The solid cross sectional area for heat conduction, AL varies 
ot 


along the path length. The area, AL , selected was the cross sectional 
i 
area between the slots normal to the flow direction, i.e., 


A, = (0.0095) (0.0022) = 0.0000209 in.? 
, | 


Just as L= a L, where L is the total matrix flow length and L. 


is the flow length between two successive perforations; so does L. = au, 
| i 
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and AL = >. A, . In evaluating Avo a conduction area reduction ratio is 


1 e 4 e 
’ cross section a of perfo ons; 
employed and is defined as £ _ £ross cti al area of perforations; 
cross sectional area , 


therefore, Aya is the solid cross sectional area multiplied by 1 -¢ 


ie.» = 
moe la e ) 
. — 0155) (.0022){ _ 2 
1.6058 C (.025) (. 0022) | ~ 0. 61024in 
= 0.00424 ft. 
160/409 
«_ -0089" T .O161" 
.0019" 


863" 


pe 





(.0063) (.025) = .0001575 in.* 


ach . | 
A_, = (.0019)(.0161) = .00003059 a 
Y = .1942 

s = .8058 
L_. = 0.3351, same as 160/40TV 
by 

S = .644 

2 

A. =A, (2-3) = -002887 ét. 


96 





125M 


ke 0079" 


—>| 
d = 00366" 
na é (5 of 





A SS 0979 sin é0° 
es WS =.006S4" 
“AQ, = (00684) (.0079) = .00005404 in.* 
A = a (.00366)- = ,00001052 in.” 
hole = 
+ “Anole = .00001052 = 0.19466 | 
A ~00005404 
ch 
s = ,8053 
le 0079" __ 
ae - —— 


Ly = .0079 sin 60° 
=.00684," 





a7 





The equivalent conduction length, L. » is computed in the case of 
2 


round perforations assuming the path to be as shown above, i.e., 27 Tie 


—— 


where th is the radial distance equal to one-half the center to center 


« 


Gistance. Thus 


Lb, = 23 .0079)= .00827" 
es 2 
and 
ro 
oe .00684 = .82709 
i —- .00827 


A is determined as in the case of the slotted perforations. 


4 = diam. of hole x plate thickness 


center to center distance x plate thickness 


= (.00366)(.0016) = .4633 
(.0079) (.0016) 
Therefore, a =A (l1- 4 ) = 1.1679in.2 (1 - .4633) 
S 
= .6268 in.* = .00435 ft.? 
E252 


be— +0079" —_s 
a = ,0029i' 
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A, = (.00684)(.0079) = .00005404 in? 


Mole" (.00291)* = 00000665 


V = 0.12307 
s = 0.8769 
Ly, = -00827" 
Li = £82709 

Lei 
SB = 0.3684 


A= AC - $y = .00640 £t° 


E 
@ 
G) 


j~— 0197" — 


d= ,.0138" 


— {+} — oe 


L, = .0197 sin 60° 


= .01706" 
a Pie 
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(.01706)(.0197) = .00033608 fre 


2 
av (0138) = .00014957 in.” 
A 


«44504 
° 55496 


02063" 


- 82697 


- 7005 


2 
A (1- 8 ) = .00243 ft. 


100 





SI 


oc 


; 0°0OS 


S°%72 


O°cT 





voiy uedo x 





[OXON pozexOfFAdG Jo satjAAodoag otAQOBWOSH 1-9 AIAVL 


6£00°0 
6£00°0 

0¢20°O 

G20 °0*€900'°0 
¢c0 °0*£900 °0 


a0ue STG 
10}uU9d9 0F 19QuUd9 


200°0 
9T00°0 
9100°0 
9100°0 
c€200°0O 


ssouyorylL 29e1Td 


€00°0 
6£00°0 
S/10°0 
810°0X6100°0 
£10°0%8000°'0 
~ (seyout) 


oZoUeTG s{oOH 


punoi 


punoz 


puno1 


po} 70[s 


pe qo [s 





sutusdg 


dz 
WSél 
90S 
607/091 


ALOV/O9T 





odd] 


101 





APPENDIX D 
SAMPLE CALCULATIONS 

The following sample calculations are based on data obtained from 
Run #1 for the modified parallel plate matrix constructed from the 160/40Q 
perforated nickel material. The basic geometric parameters are illus- 
trated in Appendix A. 

Prior to the reduction of data, several fixed parameters must be 
evaluated. 


DIMENSIONS OF FINNED SHEETS: 3.9375" x 2.00" x .0016" 
(BEFORE FORMING) 


DIMENSIONS OF SPLITTER PLATES: 3.1875" x 2.00" x 0016" 
NUMBER OF FINNED SHEETS: 102 

NUMBER OF SPLITTER PLATES: 103 

WEIGHT OF MATRIX: 311.0979gms/453.6gm/1b = 0.68584 1b. 


MATERIAL CONSTANTS: 


k = 38.7 Btu/hr - F-= ft 
S 
c, = 0.1065 Btu/lb - F 
@ = 0.321 1b/in? = 554.7 1b/£t? 


PLANE SURFACE AREA, A*, = (102)(2)(2)(3.9375 + 3.1875) = 
144 


20.1875 £t° 
HEAT TRANSFER AREA, A, = PLANE SURFACE AREA x SOLIDITY 
A = AxXs 
= (20.1875) (.8058) 


A = 16.2667 £t* 


FRONTAL AREA, A. , = 3.141" x 3.1875" = 10.022 in? 
r : 


2 
A, = 0.06953 ft 
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It has been shown in Appendix A; however, that this can be alter- 


nately expressed as: 


Inasmuch as the perforations have no effect on the hydraulic diameter 


the value of @E rather than 3 is used. 


Dy = (4) (0.8834) = | 0.002028 ft 


(1742.1) 
Now, we are in a position to reduce the data. The following presenta- 
tion. outlines the steps required; however, as far as practicable the 
actual evaluation is performed by the digital computer program set forth 


in Appendix E. 


RECORDED DATA: 


d = 2,310" T = 69.0 °F 
@] @] 

AP. = 6.69 in H,0 Patm 7 30-140 in Hg. 

yas P= 10.15 in H,0 Chart Speed = 4 in/sec. 
P = 10.65 in HO BG =d /d = 2,310 = 0.75 
8 ” 3.08 
Py = 20.80 in H0 


DETERMINATION OF MASS RATE OF FLOW: (See Appendix A) 


2 
m = 589.81 FaYdo C ; P P 
a on A oO = 2 
Y1-@ ii 
4, 
where: C=C, + AC | if from reference (15] 
N 
R(d) 
ee = CE aw - P /13.6) 0.4912; T, = t) CF) + 459.7 
Y =1.0 Fig. 40b of Reference fie] 


“f = 1.0 Fig. 30 of Reference [16] 
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For E 


0.75 and d = 3.08", C= 0.60691 and AC = 0.03839, respective- 


ly 
P. = | 30.140 - 20.80 | «x 0.4912 = 14.0537 psi 
13.6 

SOLID MATRIX CROSS SECTIONAL AREA, A., 

A, = 102 x 3.9375" x .0016" + 103 x 3.1875" x .0016" 

= 1.1679 in? 

A, = 0.00811 ft? 
FREE FLOW AREA, A> = FRONTAL AREA - SOLID CROSS SECTIONAL AREA 

A. = A; - A, 

A, = 0.06142 £t? 
CONDUCTION AREA CORRECTED FOR EFFECTS OF PERFORATIONS, AL 

A, = A (1 - 3 ) where 3 is the conduction area reduction 

? ratio (See Appendix C). 
= 0.00811 (1-0.644) 

A, = 0.00288 £t? 

POROSITY, p, = MATRIX FLOW VOID VOLUME = Vm - Vs =1 - Vs 
MATRIX VOLUME Vn Vm 
p= 1 - 0.00135 = 0.8834 
0.01159 
COMPACTNESS, @ = ‘TOTAL HEAT TRANSFER AREA 
MATRIX VOLUME 
= A - 16.2667 = |1403.75 pen eee 
Vm 0.01159 
Note that E includes the effect of area reduction due to perfora- 
tions. 


(S for an unperforated surface would therefore simply be the plane 


surface area divided by the matrix volume. 
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@=Ae = 20.1875 = 1742.1 ee7/£t? 
Vm 0.01159 


THE HYDRAULIC DIAMETER, Dy =4r, = 


4 FLOW CROSS SECTIONAL AREA 
WETTED PERIMETER 


T. = 69.0 + 459.7 = 528.7 °R 


m = (589.81 (2.310)? C 





6.69 x 14.0537 






(ae 528.7 
1 - (0.75) 

m= 1605.242 C 
assume C = 0.6000 

then M = 963.145 lb/hr 

c=C, + Ac | i 

“R(d) where NR(d,) = (4) (12) Mm 
(3.1416) (3.08) 
= 4,961 m 


pP is evaluated at temperature ty by means of the linear approximation 

Pt = 0.0395 + 0.64167 x 10° t. = 0.0395 + 0.64167 x 10°* (69.0) 
= 0.04392 lb/ft - hr 
therefore, 4 
(10 ~) (0.04392) 

C = 0.60691 + 0.03839 (4.961) (963.145) 
C = 0.6104 
This value of C is compared with the assumed value of C and successive 


iterations are performed until the two values are equal. The iterative 


process shows that a value of C = 0.6172 satisfies this requirement and 


m = 990.756 lb/hr 
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MATRIX REYNOLDS NUMBERS: 

Two values for Reynolds number are required, one for the isothermal 
flow friction factor evaluation and second for the heat transfer evalua- 
tion. The two values differ only in so far as the difference in the 
absolute viscosities. Pe is evaluated at orifice temperature (t,) 
whereas pe (heat transfer) is evaluated at average bulk fluid tempera- 


ture which is assumed to be at to + 10 °F, 





2 
Nee = 4L m = (4) 12 (990.756) = [924.54 
A pe (16. 2667) (.04392) 
No, = 4.Lm =N,. Me = 924.54 _(.04392) = 
af A Dy "t i (04456) 


CONDUCTION PARAMETER: 


r= k. Ak = (38.7) (0.00288) = 0.00282 
L wh c. (2/12) (990. 756) (. 24) 
Ax = eta ; _L_ = 0.3351 from Appendix C. 
Ly, Ly 


therefore Lo = (.00282)(.3351) = 0.000944 


MAXIMUM SLOPE: 
The slope of the cooling curve was determined by the method pres- 
cribed in Appendix A from Figure 36. Slopes were determined for each of 


the curves and the average value was used. 


MAXIMUM SLOPE = C. -2 = C.-x SLOPE x CHART 
d max — 
c 2 C SPEED 
C. =We = (0.68584 1b) (0.1065 Btu/1lb °F) = 0.07304 
Ss 
C= Mc = (990.756 1b/hr)(0.24 Btu/lb °F) _1 hr_ = 0.06605 
3600 sec 


106 





SLOPE = 0.14357 (See Figure 36) 


whence, 


MAXIMUM SLOPE = (.07304) (0.14357) (4) = 0.63507 
(.06605) 


COLBURN _} FACTOR: 
Enter Table XI with the values of A, and MAXIMUM SLOPE to deter- 
mine Neu = 4,05. Linear interpolation from this table is sufficient. 
j=Np, Ac N 2/3 from Appendix A, 
— PR 


2/3 
Nor = 0.8017 - 0.82353 x 1074 to using a linear approximation; 
therefore 
j = (4.05) (0.06142) (.7960) = |0.01217 


(16. 2667) 
FANNING FRICTION FACTOR: 
From Appendix A, the friction factor is: 
2\{ Cu 
= | 28 Qm (4E)- Kecise) = Ss (te | 
{ G* ( e e sue ( p) as 


Patm = 20-140 in Hg x 0.4912 poi = 14.8048 psi 
n ng 


P = 10.65 in H,0 x 5.204 psf = 55.42 psf = 0.3849 psi 


in H,0 

P = 10.15 in H,0 x 5.204 psf = 52.82 psf = 0.3668 psi 
m 2 fa Wo" 
P) = Pi. - Py = 14-8048 - 0.3849 = 14.4199 psi 
P, = P}- AP, = 14.4199 - 0.3668 = 14.0531 psi 
Beret oP) + PF, = 14,2365 psi 

aaa 

: 3 

Prean = (Puean) (144) = (14.2365) Gab = 0,072681b/ft 

R(t, + 459.7) (53.35) (528. 7) 
Bs) @ = 9908756 = 16130N84Ebnr - fE> —Eaeaelelby Gea — fe“ 

A 0.06142 


Cc 
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r 
H = 0.000507 = 9.00304 
i. eae 


- = 0.48 and K = -0.33 from Figure 38. 


2 
(l+p) =lt (0.8834) ° = 1.7804 


f = 2) (32.2) (0.07268) (52.82) - (0.48 - 0.33) - .3668 
(4.481) 14.2365 


(1.7804) | (0.00304) 
= (12.3126 - 0.150 - 0.04587) (0.00304) 


£ = 0.03686 


The thermal proper tes yaa i , were obtained by linear 
interpolation of data from reference [6 | - The values of kK. and C. 
were obtained directly from references | 5] and [4 | for nickel and 


stainless steel respectively. 
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APPENDIX E 

FORTRAN 60 DATA REDUCTION COMPUTER PROGRAM FOR A CDC 1604 DIGITAL COMPUTER 

A computer program was formulated in order that as much data reduction 
as practicable would be performed by the digital computer. (See Table E- 
Le 

It was hoped that this program together with a subroutine would produce 
Colburn j factor directly. An unsuccessful attempt was made to fill out 
table values of N,,, (TABLE XI) by curve fitting using Aitken's iterative 
method for fitting an nth degree Lagrange polynomial to the n + 1 given 
points - allowing for unequal spacing. As a result of reference [8 | , this 
approach was abandoned; however, the computer program presented is still 
helpful in that values of conduction parameter and maximum slope are end 
results. 

Entering TABLE XI with these parameters yields a value of Now? which 
when multiplied by Cl (Cl = N,2/3 A_/A) produces the Colburn j factor. 
A linear approximation for values of viscosity and So Me as functions 
of temperature is included. An added convenience in the computer program 
is the evaluation of the orifice constant, K. Using [15 | » an assumed 
value for C (a slowly varying function) was made. The correct value for 
C was then evaluated by an iterative process by the computer. With the 


value of C determined, K (a rapidly varying function) was found from the 


following simple relationship: K = C 
vi-e* 
where © is the ratio of the orifice diameter to the inside pipe dia- 
meter. 
With the computer program presented, an unlimited number of matrices 
can be reduced by simply changing the first card of the program. (L = 
number of matrices for which data is to be reduced). 
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The following is a list of the input data for the FORTRAN 60 program: 


(See sample input data in TABLE E - 2) 


CARD NO, 


30 


40 


10 


LG at 2 


AS 


RH 


POR 


AFL 


AHT 


DELPO 


CS 


DELH 


HS 


SLO 


TEMPO 


Matrix length, ft 

Solid conduction area, et? 

Hydraulic radius, ft 

Porosity 

Free flow area, et? 

Matrix heat transfer area, et? 

Matrix weight, lbs. 

Exit flow coefficient 

Entrance flow coefficient 

Solid thermal conductivity, Btu/hr ft °F 
Matrix material specific heat, Btu/lb °F 


Ratio of L/L, (Flow length/Equivalent conduc- 
tion length) 


Number of runs 

Diameter of orifice plate opening, inches 
Pressure drop across orifice, inches Fi 

Static pressure upstream of orifice, inches H,0 
Atmospheric pressure, inches Hg 

Ratio of orifice diameter to inside pipe diameter 
Chart speed, sec/inch 

Pressure drop across matrix, inches H,0 

Static pressure at matrix inlet, inches ae 

Slope of generalized cooling curve, (inches) 7} 


° 


Fluid temperature at orifice (t,); F 
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The output from the program is as follows: 


EMDOT 


RNUMP 


RNUMH 


CAPFL 


CAPS 


CONDPAR 


SLOMAX 


CONPARE 


G 


RNUMFR 


FFR 


ORIFICE K 


PRANDTL NO. 2/3 (Prandtl Rumen) 


Cl 


UFR 


ESTD 


STDH 


Mass rate of flow (m), lb/hr 

Reynolds number for the pipe 

Reynolds number for heat transfer 

Thermal capacitance of the fluid (C), Btu/sec °F 
Thermal capacitance of the solid (C,); Btu/ °F 
Conduction parameter ( A ) 

Maximum slope of generalized cooling curve 


Equivalent conduction parameter (A, ) 


Mass flow velocity (m/A.), lb/sec - ft? 


Reynolds number for friction 
Fanning friction factor 


Orifice constant, K 
/3 2/3 
’ Nor 
Parameter = N ee x A./A 
Pr < 
Fluid viscosity, M , lb/hr-ft 
(evaluated at fluid temperature at orifice) 


Flow friction power per unit area (E._), HP/£t2 


STD 
Parameter for evaluating heat transfer power per unit 


area (hgpp), Btu/hr ft* °F, (STDH x N,, = bgp) 
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PROGRAM SINGBLO 
L=] 
K=] 
30 READ 300sXML»2ASsRHsPOR sAFL sAHT 3 WM 
Ewe FORMAT (7F10.5) 
40 READ 4OUsEK.CKsSKsCMaXLRIN 
400 FORMAT(5F 10.5312} 
M=] 
10 READ 10U,DO,DELPO,PO;ATMP,SETA 
eee ORMAT(5F10.5) 
READ 101,CS;DELH »sHS»SLOsTEMPO 
meme ORMAT(5F10.5) 
P3=(ATMP-(PO/1346))*e4912 
TO=TEMPO+459./ 
UFR=0.0395+0 -64167E-4* TEMPO 
to ener LO. 0 
UH=0.0395+0 -64167E-4#TB 
WRITE OUTPUT TAPE 439 
9 FORMAT(1H1) 
WRITE OUTPUT TAPE 4545” 
4 FORMAT(QOH RUN NOe I1L2///) 
Miemeeteve oe (be tA. /5)--00000T)2Zis27s12 
Wemeoor (be PA—-.64)-.000001 ) 26326313 
MOMamreor (oe fae 59)—.OC O00! ) 25325514 
Ween eor (bETA—~e40)-.0U00001 ) 24324515 
ieee oor (O06 7A—2-25)—-.00000C] ) 23523516 
Memeo bor (be fA-.15)—--000001) 22922517 
meso TO 80 
21 X2=0.59225 7 
DELX=0.01916 | a 
ASSUM=0.6000 
Go TO 26 
Pex 2-0-5911 1 
DELX=0.01691 
ASSUM=0.6000 
SO 1072s 
23 XK2=0-59184 
Peex=0.01622 
ASSUM=0.6000 
“Go 10 28 us 
24 X2=0259448 
DEEX-0O.0186/7 
ASSUM=0.6000 
"Gemnro. 28 


s TABLE E-1 
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Ze 


26 


ral 


28 


eel 


29 


60 


X2=0.59850 

Welk —-0.02125 

ASSUM=0 -6000 

Goro 2&8 

R= 0.005975 

DEEX=0.02595 

ASSUM=0.60C0 

Gee1O 28 

X2=0-60691 

DELX=0.03839 

ASSUM=0.6000 

COR=ASSUM 

NC=0 

Bleie/SOORITF(1e«-BETA**4) 
NC=NC+1 

IF (NC-4)29529560 
EMBbOT=589.81*#*DO*DO*SQRTF(DELPO*P3/7T0)*COR*#B1 
RNUMO=(4.961*EMDOT/UFR)/BETA 
CORA=X24+DELX*SQRTF(1G000./RNUMO ) 
WRITE OUTPUT TAPE 4,8 5,CORA 
FORMAT (6H CORA=F10.5//) 
COR=CORA 

SO, TOna1 
PR=0.8017-.82353E-4* TEMPO 
CL=PRHAFL/AHT 

CP=0.24 

RNUMP=5 O9*# EMDOT/UFR 
RNUMH=(4-0*XML#EMDOT )/ (AHT *UH) 
RNUMFR=RNUMH*UH/UFR 
CAPFL=EMDOT#*CP/ 3600.0 “= 
CAPS=WM*CM . 

eIperiAx=( CAPS*#SLO) /(CAPFL*CS) 
GONDPAR=(SK*#AS)/(XML*EMDOT*CP) 
CONPARE=CONDPAR*¥XLR 

B=CK+EK 
C=1.20+POR**2 
G=EMDOT/(AFL*3600.0) 
Peee-DELH*.03613 

BeeiioteO3613 ° 

PA=ATMP*.4912 : 

PLT=PA—-PS 

P2=P1-DELP 

Pret P1+P2)/2. 
RHOM=(PM¥144.0)/(53- 3%(TEMPOt460.0) ) 


TABLE E-1 
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FF 1=((64.4*RHOMEDELP#¥144.0)/(G#G)) 
FF2=C*#*DELP/PY 
Pe S=FPPi=Fr2-B 
Pate 3* RAZ XML 
CON T=CCOR3 1 
STDH=(0202195*RNUMH*¥C1)/(4.*RH) 
eo ptaete LlE—/*( 1. 5 / Ws #RH) VERE eRe TRUER Cmte 
WRITE OUTPUT TAPE 43575,STDH;ESTD 
7 FORMAT( 19H CONSTANT FOR HSTD=F102553X95H ESTD= Ore 5/11) 
WRITE OUTPUT TAPE 451 ,sEMDOT »sRNUMPsRNUMHs CAPFL ,CAPS 
1 FORMAT( 7H EMDOT=F10.533Xs7H RNUMP=E20¢853X9s7H RNUMH=F10.593Xs5 
vinee@re rm le—— 10. 553Xs6H CAPS=F10.5///7) 
WRITE OUTPUT TAPE 432sCONDPAR »SLOMAX , CONPARE 
2° FORMAT(QOH CONDPAR=F102593Xs8H SLOMAX= ee: 593X%99H CONPARE= POs 57 7a 
WRITE OUTPUT TAPE 45935,GsRNUMFRSFFRCONS 
BORO RMAIt sd G=—r-1l0.553X+8H RNUMFR=FIC. ere SH FFR=P 10455 3Xs1 1 Ones 
Wesek—=— 10. 2// / } | 
WRITE OUTPUT TAPE 4555sPRsC1SUFR 
5 FORMAT( 16H PRANDTL NOe2/3=F106 5 aeea Cl=F L065 3X55 VER=f eee 
M=M+] 
IF (N+M) 50310310 
50 K=K+1 
Tet G=K )205320:530 
80 PRINT 81 
Sheena ( 20ReORTe ICE CONSTANT ERROR EXISTS) 
GO 10620 
AO 5) OF 
END 
END ‘ 


TABLE E-1 
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SAMPLE INPUT DATA FOR PROGRAM SINGBLO 


0.1667 
ma ) 32 
Ze 310 
025 
22310 


2.310 
025 | 
1.971 
250 
1.540 
250 
1.540 
050 
1.232 
050 
0.775 
1.00 
00775 
1.00 
0.775 
2200 
0.462 
4.00 


e00640 
051 
6278 
9.30 
4.40 
7e15 
2074 
5215 
382 
3240 
6289 
2233 
2091 
1242 
2e/5 
0-80 
6200 a: 
0241 
5.98 
0-240 
1.78 
0622 
2260 
02080 


e00049 
38e7 
2145 
l2e2 
15245 
8 e 30 
5280 
7e¢10 
3250 
450 
2206 
2260 
0095 
1240 
0243 
0268 
0218 ~ 
0267 
0-178 
0235 
02075 
0.15 
02025 


TABLE 


©8542 

©1065 
302140 
© 17768 
302120 
015259 
30.2120 
012882 
302120 
e21040 
302120 
017379 
302110 
e 14039 
30.2110 
009232 
302110 
e10170 
302110 
© 10316 
302110 
©1032] 


302110 


007513 


05927 
e82709 


075 
685 
075 
7120 
075 
730 
064. 
7220 
250 
7207 
250 
7207 


"240 


7245 
025 
69.0 
025 
69.0 
025 
7020 
e15 
7020 


17270242 
1p Nis 


04844 
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